Stimulated Raman adiabatic passage analogs in classical physics 
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Stimulated Raman adiabatic passage (STIRAP) is a well established technique for producing 
coherent population transfer in a three-state quantum system. We here exploit the resemblance 
between the Schrodinger equation for such a quantum system and the Newton equation of motion 
for a classical system undergoing torque to discuss several classical analogs of STIRAP, notably the 
motion of a moving charged particle subject to the Lorentz force of a quasistatic magnetic field, 
the orientation of a magnetic moment in a slowly varying magnetic field, the Coriolis effect and the 
inertial frame dragging effect. Like STIRAP, those phenomena occur for counterintuitively ordered 
field pulses and are robustly insensitive to small changes in the interaction properties. 



PACS numbers; 32.80.Qk,32.80.Xx, 33.80.Be, 91.12.Hg 



I. INTRODUCTION 

The production of complete population transfer be- 
tween the quantum states of a three-state chain by means 
of stimulated Raman adiabatic passage (STIRAP) has 
become one of the staples of contemporary quantum- 
state manipulation [1-4]. In this technique, a pair of 
temporally delayed laser pulses P and S, of frequencies 
ujp and Lus respectively, act to alter the state vector ^'(t) 
in a three-dimensional Hilbert space from initial align- 
ment virith state ipi to final alignment with state tp^ . The 
P field links state tpi, of energy Ei, to the intermediate 
state 1^2, of energy E2, and the 5' field links this to the 
final state ips, of energy E3. By applying the S field prior 
to the P field (so-called counterintuitive ordering), and 
maintaining two-photon resonance along with adiabatic 
evolution, the population transfer occurs without state 
■02 acquiring, even temporarily, any population. 

In many classical systems the equations of motion com- 
prise three coupled equations which can be cast into the 
form of a torque equation, i.e. an equation of motion in 
which the force acting on a vector is always at right an- 
gles to the vector. The behavior of a gyroscope acted on 
by gravity is a familiar example. Such a torque equation 
occurs in quantum optics as the well-known Bloch- vector 
representation of the behavior of a coherently driven two- 
state quantum system [5, 6]. It is less well known that 
a torque equation applies to the three-state system of 
STIRAP [7]. 

In this paper we exploit the occurrence of a torque 
equation in these several different areas of physics — 
quantum mechanics, classical mechanics, classical elec- 
trodynamics and general relativity — to discuss ways in 
which classical motion can be altered adiabatically us- 
ing two sequential but overlapping pulsed interactions. 
In particular, we will note the analogy with the three 
Cartesian coordinates of a moving charge in electric and 
magnetic fields, as described by the Lorentz force [8], 
and the orientation coordinates of a magnetic moment 



in a magnetic field, as described by the Landau-Lifshitz- 
Gilbert equation [9, 10]. 

We note that several authors have discussed analogies 
between three-state quantum systems and classical sys- 
tems. These similarities include an analog with the mo- 
tion of a classical pendulum [11] and an analog with clec- 
tromagnetically induced transparency [12]. 

This paper is organized as follows. In Sec. II we 
present the basic mathematics of STIRAP, and cast the 
equations into the form of a torque equation. In Sec. Ill 
we present the equations of motion for a charged particle 
subject to a Lorentz force, with specialization that makes 
these identical to the STIRAP torque equation. Specif- 
ically, we consider a particle that initially moves in the 
z direction, acted upon by a sequence of two magnetic- 
field pulses. The first of these has only a z component 
(thereby producing no change in the motion), while the 
second has only an x component. The resulting particle 
motion is in the x direction; there is never any component 
in the y direction, despite the presence of the x-directed 
magnetic field. In Sec. IV we discuss three more ex- 
amples of classical physics where STIRAP-like processes 
can occur. These include the reorientation of magnetiza- 
tion, the Coriolis effect and the general relativity effect of 
inertial frame dragging. Section V presents a summary. 



II. STIRAP 

The basic equation of motion governing STIRAP is the 
time-dependent Schrodinger equation, in the rotating- 
wave approximation (RWA). Expressed in vector form, 
this reads 



ih—c{t) = H(i)c(i), 



(1) 



where c(t) is a column vector of probability amplitudes 
Cn{t) and H(t) is the RWA Hamiltonian matrix [13-15]. 
In the example of STIRAP there are three basic quantum 
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states — ■01, ^2 and ips — linked as a chain 1 
and H(i) is a 3 X 3 matrix, 



rjp(t) 
fip(t) risit) 

llsW 



(2) 



Here the two slowly varying Rabi frequencies flp{t) and 
ils{t) parameterize the strengths of the pulsed P- and 
S'-field interactions; they are proportional to dipole tran- 
sition moments dij and to electric-field amplitudes (t) , 



-driSpit), - ^d23£s{t), (3) 



and hence they vary as the square root of pulse in- 
tensities. We take these to be real-valued functions of 
time. We have here assumed not only two-photon reso- 
nance h\ujp — ujs\ = \E3 — Ei\, as required for STIRAP, 
but also single-photon resonances: \E2 — Ei\ = hiup, 
\E2 — E^l = huJs- 

The quantum evolution associated with STIRAP is 
most easily understood with the use of adiabatic states, 
i.e. three instantaneous eigenstates of the RWA Hamilto- 
nian. For the assumed condition of two-photon resonance 
between ipi and V's i one of these adiabatic states has no 
component of state ?/'2; it is a dark state. In the absence 
of the P field this adiabatic state coincides with state ipi , 
while in the absence of the S field it is aligned with state 
■03. By maintaining adiabatic conditions the state vec- 
tor ^{t) follows this Hilbert-space motion, thereby pro- 
ducing complete population transfer from state -01 to '03. 
Moreover, the motion of the state vector remains entirely 
in a two-dimensional subspace of the three-dimensional 
Hilbcrt space. This restriction allows a simple descrip- 
tion of the dynamics as a torque equation. 

The instantaneous eigenvectors ^k{t) of the matrix 
H{t), defined as H(t)$fe(i) = fefc(i)$fe(t), are 



V'2 + 03Cosi9(t)], (4a) 
ipiCOSiSit) -ip3sm-d(t), (4b) 
^ [i'l sin^it) - VJ2 + 03 cosi?(<)] . (4c) 



Here the time-dependent mixing angle d{t) is defined as 
the ratio of interaction strengths. 



npjt) 



(5) 



The adiabatic state ^o{t) is particularly noteworthy: it 
has a null eigenvalue and it has no component of state 
02 — it therefore does not lead to fluorescence from that 
state; it is a dark state [1-4, 16-19]. The construction 
of all three adiabatic states varies as the pulse sequence 
alters the mixing angle, but ^o{t) remains at all times 
within a two-dimensional Hilbert subspace. It is this 
property that we exploit for our classical analogies. 

The two alternative pulse orderings, S ~ P and P ~ S, 
lead to qualitatively different results. 



A. Counterintuitive pulse order: STIRAP 

The STIRAP mechanism relies on maintaining a con- 
tinuing alignment of the state vector \l/(t) with the dark 
state ^o{t), and having this adiabatic state initially 
aligned with state ■01- To have this initial alignment it 
is necessary that the S field act first. Because this field 
has no interaction linkage with the initially populated 
state 01 it does not directly produce population transfer. 
Thus the 5'-before-P sequence is termed counterintuitive. 
With this ordering the mixing angle ■d{t) and the adia- 
batic state $o(i) have the behavior 



01 <f>0(i) -^3- 



(6a) 
(6b) 



That is, the S — P pulse sequence rotates the adiabatic 
state ^a{t) from alignment with the initial state 0i to 
alignment with the target state "03. If the motion is 
adiabatic, then the state vector ^{t) follows this same 
Hilbert-space rotation. The result is complete popu- 
lation transfer. The condition for adiabatic evolution 
amounts to the requirement of large temporal pulse ar- 
eas, A„ = /_~^ Qkit) dt {k = P, S) [2-4], 



Ap > 1, > 1- 



B. Intuitive pulse order: oscillations 



(7) 



The intuitive sequence P—S produces populations that 
display Rabi- like oscillations [20] . This behavior is read- 
ily understood by viewing the construction of the adi- 
abatic states for very early and very late times. Let 
the state vector coincide with state 0i at time t — > — cx), 
when the S field is absent. Then this initial state has the 
construction [cf. Eqs. (4)] 



*(-oo) = 01 = ^[$+(-CX)) + $_(-CX))]. 



(8) 



Later, towards times t — > -|-c» when only the S field is 
present, this construction becomes 



o) = ^[e-''^+(°°)$+(-oo)-f e-'"^~(°°)$_(-cx))], 

(9) 

where 0fe(i) = /gefc(t')di' {k = +,— ). The oscillatory 
phases lead to oscillations of the populations [20] , 

Pi =0, P2 = sin^ iA, P3 = cos^ lA, (10) 

where A is the rms pulse area, 

A = 0+(cx)) - 0_(oo) = f ^Jnj,{t) + nl{t) dt. (11) 

Thus, only for certain values of the pulse area (general- 
ized TT-pulscs), it is possible to obtain complete popula- 
tion transfer from state 0i to state 03 with a resonant 
P ~ S pulse sequence. 
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C. The STIRAP torque equation 

The three-state Schrodinger equation driven by the 
Hamiltonian (2) has, with a redefinition of variables, 



the form 



_d 



-C3(t), R2{t) 



-iC2{t), i?3W=Ci(t), 

(12) 



Rlit) 




R2{t) 























Ri{t) 

R2{t) 

Riit) 



(13) 

The symmetry of these coupled equations allows us to 
write this equation of motion as a torque equation [7] , 



—R{t) = n{t)xR{t), 
by introducing the angular velocity vector 



(14) 



(15) 



Like torque equations in classical dynamics [21], this 
equation says that changes in a vector, R(t) in this case, 
are produced by a force J7(t)xR(<) that is perpendicular 
to the vector. Such an equation occurs in the description 
of two-state excitation, where it gcometrizcs the Bloch 
equation [5]. In all cases it describes an instantaneous 
rotation of the vector R(t), in the plane orthogonal to 
the direction of the angular velocity vector at an instan- 
taneous rate |fl(t)| = y/^pit)^ + ^sit)^- 

The dark state (4b) written for the vector H{t) gives 
the dark superposition 



Roit) 



np{t)Riit) + nsit)R3it) 

m)\ 



(16) 



The STIRAP evolution, regarded as the solution to 
the torque equation (14), comprises the following mo- 
tion. Initially R(i) points along the z-axis and n{t) 
points along this same axis. Because these two vectors 
are collinear the vector R(i) does not move. The intro- 
duction of the P field rotates fl{t) in the xz-p\a,ne toward 
the X axis. Because this motion is adiabatic, it causes the 
vector R(t) to follow. In the end, both R(i) and ^{t), 
remaining collinear, are aligned along the x-axis. 

We note that the conservation of the length of the vec- 
tor R(i), which follows from the torque equation (14), 
is equivalent to the conservation of probability ensuing 
from Eq. (1). 

The following sections will note several classical sys- 
tems that are governed by a torque equation, and will 
discuss the analogs of STIRAP motion for these systems. 



III. STIRAP IN LORENTZ FORCE 

A charged particle moving in a magnetic field is altered 
by a force that is perpendicular to both the velocity and 




FIG. 1: (Color online) Evolution of the magnetic vector (top 
frame) and the velocity vector of the moving charged parti- 
cle (bottom frame). For counterintuitive order of pulses (top 
frame) the initial velocity changes from alignment with z axis 
to alignment with x axis without having any components 
on the y axis. We have assumed Gaussian magnetic pulse 
shapes, B4t) = Bo e-^'+^^^^' B,(t) = Bo e'^'-^^^^" 
with Bo = 20T~^ and r = -1.2T. 



the magnetic field. Let the particle have a mass m, a 
charge q and a velocity v(t) = [vx{t),Vy{t),Vz{t)]'^ . Let 
the magnetic field be restricted to components in the 
x2-plane, B{t) = [-B^{t),0, B^it)]'^ . The Lorentz force 
acting on the particle is x B. The Newton's equation 
of motion therefore appears as a torque equation for the 
particle velocity, 



m — V = — oBxv. 

dt 



(17) 



Drawing an analogy to STIRAP, we write down a dark- 
velocity superposition Vo{t) of the velocity components 
Vx{t) and Vz(t), 



Voit) 



B,{t)v,{t) + B,{t)v,it) 
|B(t)| 



(18) 



When Bz (t) precedes B^ (t) then the dark velocity super- 
position Vo(t) has the asymptotic values 



w^(oo). 



(19) 



Thus if initially the particle travels along the z direction, 
v{t) = [0,0, w]"^, we can direct the velocity into the x 
direction by applying first a field in the z direction and 
then slowly rotating this into the x direction as shown 
in Fig. 1. The initial field, being in the direction of 
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and the length L, 



Pulse Delay (in units T) 

FIG. 2: (Color online) Numerically calculated velocity com- 
ponents Vx{t) (dashed black line), Vy{t) (doted red line) 
and Vz{t) (solid blue line) as a function of the pulse delay. 
We have assumed Gaussian magnetic pulse shapes, B^it) — 
Boe-(*+"/2)'/^', B^t) = Boe-"-^''2>'/^' with Bo = 40r-\ 



motion, has no effect; in this sense the pulse sequence is 
counterintuitive. 

If the initial magnetic field is in the x direction (intu- 
itive pulse order), then the charged particle, which travels 
initially along the z axis, will be subjected to a Lorcntz 
force and will begin a Larmor precession in the j/z-planc. 
Then, as the B-field switches from x to z direction, the 
particle precession will turn into the xy-plane. The fi- 
nal velocity of the particle depends on the value of the 
accumulated precession angle A [cf. Eqs. (10) and (12)]: 
v(+oo) = [cos2(A/2),sin2(A/2),0]^. 

These features are demonstrated in Fig. 2, where the 
velocity components of the charged particle are plotted 
versus the delay between the magnetic pulses Bx{t) and 
Bz (t) . A fiat plateau of high values of Vx is observed for 
negative delays (counterintuitive pulse order), whereas 
oscillations between Vx and Vy occur for positive delays. 
Note that the final value of Vz does not depend on the 
sign of the delay r [22] . 

The equation of motion (17) holds only for quasistatic 
fields. To implement the desired STIRAP analogy we can 
use a spatial arrangement of the magnetic field such that 
the components appear to the moving particle as two 
sequential, but overlapping magnetic fields. This spatial 
geometry, viewed in the reference frame of the particle 
allows us to write the magnetic field as time-dependent 
without any associated electric field. 

Following quantum-optical STIRAP, we conclude that 
the condition for adiabatic evolution in this Lorenz- 
STIRAP is a large value of the accumulated precession 
angle A (which is the pulse area in quantum-optical STI- 
RAP). Within the above mentioned spatial arrangement 
(with a characteristic length L), the adiabatic evolution 
condition sets an upper limit on the charged particle ve- 
locity V, or lower limits on the peak magnetic field Bq 



mv ^ qBoL. 



(20) 



IV. OTHER EXAMPLES 
A. Magnetization 

A magnetic field H acts to turn a magnetic moment 
M{t) = [Mx{t),My{t),Mz{t)]'^, with a force that is al- 
ways perpendicular to M(i). The system dynamics is 
expressible again as a torque equation, 



dt 



M{t) = 7M(<) X H(<), 



(21) 



where 7 is the gyromagnetic ratio. This is the homo- 
geneous Bloch equation for magnetization with infinite 
relaxation times [6]. It is also known in the literature as 
the undamped case of the Landau-Lifshitz-Gilbert equa- 
tion [9, 10]. The dark superposition for the magnetic 
moment reads 



Mn 



Hx{t)Mx{t) + Hz{t)Mz{t) 
\M{t)\ 



(22) 



When H(f) = [—Hx{t),0, Hz{t)], and the magnetic com- 
ponent Hz{t) precedes the magnetic component Hxit), 
the dark magnetic moment M^it) has the asymptotics 



A/.(-«)) 



Mx{+oo). (23) 



magnetic moment 
direction, M(t) = [0,0, A/]^, we can 



Thus if we start with the initial 
pointed in the z 
change the direction of the magnetization from the z axis 
to the X axis by applying first a magnetic pulse Hz{t) 
and then a magnetic pulse Hx{t) (counterintuitive or- 
der), while maintaining adiabatic evolution. Because the 
adiabatic passage is robust, this procedure is robust: it 
depends only weakly on the overlap of the two magnetic 
components and the peak values of Hx{t) and Hz{t). 



B. Coriolis effect 

In classical mechanics, the Coriolis effect is an appar- 
ent deflection of a moving object when it is viewed from 
a rotating reference frame. The vector formula for the 
magnitude and direction of the Coriolis acceleration is 



dt 



V = 2v X cj. 



(24) 



where v(t) — [vx{t)TVy{t),Vz{t)] is the velocity of 
the particle in the rotating system and u}(t) = 
[uJx{t),LOy{t),uUz{t)\ is the angular velocity vector of the 
rotating frame. This equation has the same vector form 
as Eq. (14) and therefore, a STIRAP-like process may be 
demonstrated if the angular velocity vector of the rotat- 
ing frame changes appropriately. 
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C. General relativity 

An equation of the form (14) emerges in the descrip- 
tion of the effect of general relativistic gravitational frame 
dragging, e.g. when a massive spinning neutral particle 
is placed at the center of a unidirectional ring laser [23]. 
Then the linearized Einstein field equations in the weak- 
field and slow-motion approximation lead to an equation 
for the spin of the same form as Eq. (14). 

V. CONCLUSIONS 

We have presented several examples of well-known dy- 
namical problems in classical physics, which demonstrate 
that the elegant and powerful technique of STIRAP in 
quantum optics is not restricted to quantum systems. 
The application of STIRAP to these problems, which 
appears experimentally easily feasible, is intriguing and 
offers a potentially useful and efficient control technique 
for classical dynamics. 

The first factor that enables this analogy is the equiv- 
alence of the Schrodinger equation for a fully reso- 
nant three-state quantum system, wherein the quantum- 
optical STIRAP operates, to the optical Bloch equation 
for a two-state quantum system. The second factor is the 
Feynmann-Vernon-Hellwarth vector form of the Bloch 
equation, which has the form of a torque equation, i.e. 



the force on a vector is perpendicular to the vector. 

In the Lorentz force case, the variables for the STIRAP 
analogy are velocity components. The STIRAP proce- 
dure changes the direction of the velocity from the z axis 
to the X axis with never a component along the y axis. 
The procedure has the same efficiency and robustness as 
STIRAP. The described technique for a Lorentz force is 
not only a curious and intriguing example of the adia- 
batic passage, but it also has the potential to be a use- 
ful, efficient and robust technique for magnetic shielding, 
magnetic lenses, or speed selection of charged particles. 

Applied to the equation of motion of a magnetic mo- 
ment in a magnetic field, the analogy of STIRAP offers a 
robust mechanism for changing the orientation of a mag- 
netic moment. STIRAP-like processes can also be de- 
signed in other intriguing physical situations, such as the 
Coriolis effect and the general relativity effect of gravita- 
tional frame dragging. 
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